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ABSTRACT 

We use the best available X-ray data from the intermediate polar EX Hydrae to study the cooling-flow model often applied to interpret 
the X-ray spectra of these accreting magnetic white dwarf binaries. Eirst, we resolve a long-standing discrepancy between the X-ray 
and optical determinations of the mass of the white dwarf in EX Hya by applying new models of the inner disk truncation radius. 
Our fits to the X-ray spectrum now agree with the white dwarf mass of 0.79 Mg determined using dynamical methods through 
spectroscopic observations of the secondary. We use a simple isobaric cooling flow model to derive the emission line fluxes, emission 
measure distribution, and H-like to He-like line ratios for comparison with the 496 ks Chandra High Energy Transmission Grating 
observation of EX Hydrae. We find that the H/He ratios are not well reproduced by this simple isobaric cooling flow model and 
show that while H-like line fluxes can be accurately predicted, fluxes of lower-Z He-like lines are significantly underestimated. This 
discrepancy suggests that an extra heating mechanism plays an important role at the base of the accretion column, where cooler ions 
form. We thus explored more complex cooling models, including the change of gravitational potential with height in the accretion 
column and a magnetic dipole geometry. None of these modifications to the standard cooling flow model are able to reproduce the 
observed line ratios. While a cooling flow model with subsolar (0.1 O) abundances is able to reproduce the line ratios by reducing 
the cooling rate at temperatures lower than 10^'^ K, the predicted line-to-continuum ratios are much lower than observed. We 
discuss and discard mechanisms, such as photoionization, departures from constant pressure, resonant scattering, different electron- 
ion temperatures, and Compton cooling. Thermal conduction transfers energy from the region above lO’ K, where the H-like lines 
are mostly formed, to the cooler regions where the He-like ions of the lower-Z elements are formed, hence in principle it could 
help resolve the problem. However, simple models indicate that the energy is deposited below 10* K, which is too cool to increase 
the emission of the He-like lines we observe. We conclude that some other effect, such as thermally unstable cooling, modifies the 
temperature distribution. 
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1. Introduction 

Cooling flows have been postulated to occur in a diversity of 
settings, from clusters of galaxies, where the low-density, intr- 
aclus ter medium loses its gravitational energy by ra diating X- 
rays (iFabian & Nulseniri97^ : |Binnev & Co^^ 198 II) to accret¬ 
ing magnetic and nonmagnetic white dwarfs (WDs) in binary 
systems. In nonmagnetic WDs.^^ such a s dw arf novae (Patterson 
& Raymond 19851) and symbiotics (ILuna & Sok ol oskil 20071) a 
cooling flow is present in the boundary layer of the accretion 
disk, where material has to dissipate its rotational energy be¬ 
fore settling onto the WD surface. If the WD magnetic field is 
strong enough to disrupt the accretion disk, the material flowing 
through the magnetically-dominated accretion column reaches 
highly supersonic velocities and undergoes a strong shock near 
the WD surface, heating the plasma to keV temperatures; the 
material in the postshock region cools until it settles onto the 
WD surface. 

In intermediate polars (IPs), the WD has a magnetic field 
stro ng enough to disrupt the accretion disk at t he Alfven radius 
Rin (iRevnivtsev et alJl201 ll : ISemena et^ 20141) where mag¬ 
netic pressure balances the ram pressure of the accretion flow. 
The WD’s magnetic field then channels the material from the in¬ 


ner disk toward the magnetic poles. Most of the cooling radiation 
from the accretion column is emitted at X-ray wavelengths, pre¬ 
dominantly as bremsstrahlung and line emission, i.e., a cooling 
flow. 

For galaxy clusters, spatially resolved X-ray imaging obser¬ 
vations are available, and surface brightness and temperature 
profiles have been compared with the predictions of the cool¬ 
ing flow model. The observations do not show the large amounts 
of cool X-ray emitting plasma at the center of the clusters pre¬ 
dicted by the cooling flow model. This discrepancy i s known 
as the cooling flow problem dPeterson & FabianI 2006) As the 
X-ray emitting region of the accreting WD cannot be spatially 
resolved, an analogous cooling flow problem has been more 
difficult to discern. Hints of problems with the simple cooling 
flow model came from residuals in the fits of low- and medium- 
resolution X-ray spectra. However, these residuals can be re¬ 
duced by invoking complex absorption (e.g.. lRamsav et al.l2008l : 
lYuasa et alJ 201Q) Further progress has been stymied by the low 
resolution of CCD spectra and t he limited nu mber and quality of 
high-resolution grating spectra. iMukai et ^ (20031) modeled the 
high-resolution Chandra X-ray spectra of EX Hydrae (EX Hya), 
two dwarf novae, and an old nova using a solar abundance cool¬ 
ing flow model, matching the line fluxes to within a factor of two. 
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On the other hand, hts with the cooling flow model to the high- 
resolution X-ray spectra of the dwarf nova WX HyiJPerimet^ 


2003h and the likely magnetic system V426 Oph (iHomer et al 


2004) showed good agreement with the high temperature part of 
the spectrum, but underpredicted emission from the cooler part 
of the flow. This discrepancy appears to be in opposite direction 
from that seen in galaxy clusters. It is possible that there is ad¬ 
ditional heat as the gas cools to 10^ K in both cases, but that in 
EX Hya this heating simply increases the emission measure of 
the million degree gas, whereas in clusters the heating is strong 
enough to prevent the gas from cooling to lower temperatures. 

Depending on the physical properties of the binary system, 
additional heating and cooling mechanisms can play important 
roles in the energy equation of the cooling flow; e.g., gravity is a 
source of extra thermal energy when the shock height is a signifi¬ 
cant fraction of the WD radius. The cooling mechanism has been 
the subjectofman^theoreticalstudies_(ejgjAizu|l97^Cropper 
et al. ll999l : ICanalle et al.l2005F Havashi & lshidal 2014a[), and the 
derived spectral models have shown reasonable success in tests 
on low- and medium-resolution X-ray spectra. Given that the 
shock temperature is proportional to the strength of the gravita¬ 
tional potential well, its measu rement can be used to determine 
the mass of the accreting WD. Ilshida et’aD (199 ID modeled the 
X-ray spectra of IPs observed with GINGA as bremsstrahlung 
emission and used the deriv ed temperature s to determine the 
masses of the accreting WDs. ICropper et alJ ( i999D also derived 
WD masses from GINGA data, but considered additional terms 
in the energy equation of the cooling flow, such as gravity, as 
well as modifications of the observed spectra by reflection by the 
WD surface and absorption due to a warm absorber above the 
emissign,TggionJSuleimano^,e^.alJ_i2QQ5LjQiBrunschweiger 
et al. (2009) modeled RXTE an d Swif i/BAir s pectra of a sample 
of IPs using the ICropper et alJ ( 1999D model. lHavashi & Ishidal 


(2014aD combined the magnetic held geometry, varying gravi¬ 
tational potential, nonequipartition between ions and electrons, 
and nonequilibrium ionization to derive a model for the post¬ 
shock region of the accreti on column, which was rec ently tested 
on Suzaku data of EX Hya (iHavashi & Ishidal 2014b) 

Given that the cooling plasma also produces line emission, 
lEuiimoto & Ishid^ (1997D realized that the then-higher spectral 
resolution of ASCA CCD spectra could be used to determine the 
shock temperature, and hence the WD mass, through measure¬ 
ments of the line intensity ratios of the H-like to He-like ions. 
They assum ed the temperature and density structures given by 
lAizul (1973D, where cooling is due only to bremsstrahlung radia¬ 
tive losses. Although line fluxes and equivalent widths from Ee, 
Ar, S, Si, and Mg were measured. Mg lines were not used to de¬ 
termine the WD mass because they are blended with Ee L-shell 
lines at the A5CA/SIS spectral resolution. 

Are the more recent, high-quality spectra obtained with the 
Chandra High-Energy Transmission Grating (HETG) compat¬ 
ible with the model s proposed for the cooling of the accretion 
column in EX Hva? lLuna et al.l (2010) discovered high-velocity 
wings on several bright emission lines in the 496 ks HETG spec¬ 
trum of EX Hya, which they attributed to photoionized gas in the 
accretion column above the shock. Since only a small fraction of 
the X-ray spectrum originates in the photoionized preshock ma¬ 
terial, with the rest attributed to collisionally ionized plasma, a 
cooling flow model of the postshock emission should give a good 
first-order match to the observed spectrum. 

In Section|2]we describe the properties of EX Hya. We solve 
the long-standing controversy about the WD mass of EX Hya 
derived using X-ray spectral fits compared with optical/IR radial 
velocity measurements. We use the 496 ks Chandra HETG spec¬ 


trum described in Section [3] to test models of the cooling in the 
postshock region (Section|4]i. We first determine the shock tem¬ 
perature, which matches the value implied by a WD mass of 0.79 
Mq, from the least-absorbed spin-phase in Section 144] We then 
utilize the H-like to He-like line ratios as a proxy for the temper¬ 
ature and density structure of the column. We start by computing 
a simple isobaric cooling flow spectrum due to bremsstrahlung 
and line emission. We compute the emission measure distribu¬ 
tion (EMD), line fluxes, and line ratios (Section [4. 2l i. We And 
that simple isobaric cooling flow models are not able to repro¬ 
duce the observed spectrum and line ratios. Specifically, we And 
that the H-like line fluxes are well reproduced, whereas the same 
model underestimates the He-like line fluxes. This implies that 
either heat is not efficiently removed in the upper portion of the 
accretion column, or there is an extra source of heat in the lower 
portion of the accretion column. We then consider the effects 
on the EMD of a different geometrical configuration by calcu¬ 
lating a magnetic dipole model (Section [4.31) . We And that the 
line ratios predictions are mostly unaffected by the accretion col¬ 
umn geometry. In Section l4~4l we compute subsolar abundance 
models. While models with very low abundances are able to re¬ 
produce the observed line ratios, they greatly underestimate the 
observed line-to-continuum ratios. 

Motivated by recent studies of the importance of thermal 
conduction in the co oling of the intracluster medium of clus¬ 
ters o f galaxies (e.g.. lKim & Nara Yan|200 3l;^Ruszkowski & (5hl 
201 ID and WD boundary layers (TLiu et al. 2008D, we consider 
thermal conduction in the energy equation of the flow in the ac¬ 
cretion column of EX Hya (Section l43]) . While this changes the 
EMD in the correct sense, the change is too small to alleviate the 
discrepancy. We conclude that another effect, such as thermally 
unstable cooling, must be responsible for the excess emission in 
the 10^ to 10^ K range. In Section|5]we discuss our findings and 
present concluding remarks. In the Appendix we discuss various 
approximations in the simple cooling model. 


2. EX Hya 

EX Hya is an IP, a magnetic cataclysmic variable (CV) whose 
magnetic field is strong enough to channel the accretion flow 
onto the magnetic poles, but not strong enough to entirely dis¬ 
rupt the accretion disk. The orbital and WD spin periods are 
98.26 and 67.03 minutes, respectively. EX Hya has been ex- 
tensiv ely studied at optic a l, UV, EUV, and X-ray wavelengths 
(e.g., Rosen et alJ 19881: Mauchelll999l : iMukai et alJ 120^031 : 
iBrunschweiger et al.l 2009) The luminosity and accretion rate 
are estimated to be (2.6 + 0.6) x 10^^ erg s“' and (6.2 + 
1.5) X IQ-'^ (Mwn/0.5 Mo yr^' at a distance of 64.5 

pc (iBeuermann et al.l 2003D There has been a long-standing 
discrepancy between the WD mass of 0.790 + 0.026 Mq de- 
rived from the radial velocitie s of optical and IR emission lines 
(iBeuermann & Reinschl 2008D and masses of 0.42-0.66 Mq de¬ 
rived by equating the X-ray temperature to the shock tempera - 
ture expected for free fall onto the WD (e.g.. lYuasa et al.l 2010) 
Various models have been used to fit the X-ray spectrum of 
EX Hy a, and differ ent values for the WD masses have been ob¬ 
tained. [^asa^eT^ (2010) modeled the Suzaku X-ray spectrum 
including gravi ty in the flow equ ations and obtained a mass of 
0.42+0.02 M,T^: ISuleimanov et al.l (2005D modeled the RXTE X- 
ray spectrum also including gravi ty in the flow equatio ns and 
obtained a mass of 0.50+0.05 M,T^: IBrunschweiger et al.l (2009) 
obtained a mass of 0.66+0.17 Mq fitting the 5w/ft/BAT spectrum 
and using_Agjgmg^ 2 dsLi 2 Jl^l£iniiI 12 L£L 2 lii 2222 LJ^^y^^hi 
& Ishida (20141i) combined the magnetic field geometry, vary- 
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ing gravitational potential, nonequipartition between ions and 
electrons, and ionization nonequilibrium to develop a model that 
was applied t o Suzaku data of EX Hya , and obtained a mass of 
0.63^g MrT'. lFuiimoto & Ishidal l 19971) measured the ratio of the 
intensities of H-like to He-like emission lines in their study of 
the ASCA spectrum of EX Hya and obtained a mass of 
Mq. Radial velocity measurements of the X-ray emission lines 
observed in a short (50 ks) Chandra observation, combined with 
velocjtijmElitude^^_of_the_secondari^eriYedi^Vande Putte 
et al. (20031) led iHoogerwerf et al.l (2004) to derive a mass of 
0.49+0.13 Mq. Dynamical methods through optical and near- 
IR spectroscopic observations of the secondary bv Beuermann 
& Reinsch (20081), significantly improved the measurements of 
K 2 and concluded that the mass of the WD is 0.790+0.026 M©. 
The WD mass d erived by Hoogerwerf et akl (2 004) matches the 
yalue deriyed by lBeuermann & Reinschl ( 20081) using their more 
accurate measurement of K 2 . 

All the preyious determinations of the WD mass based on the 
shock temperature of EX Hya, assumed that the material in the 
accretion column was at free-fall speed from infinity. Howeyer, 
in systems such as EX Hya, the in ner radius of the accre tion 
disk can be close to the WD surface dReyniytsey et alJ 201 If) and 
therefore a factor that takes this into account has to be introduced 
in the free-fall speed equation as 


Vff - 2GMwd(— - : -— I, 

\ \RWD+/Jsh Rin/ 


( 1 ) 


where Rwd is the radius of the WD and hgt is the height of 
the accretion column. If we use Mwn -0-79 Mq, Rwn=7xl0^ 
cm, tish^O.lxRwn, and Rin=2.7xRwD dReyniytsey et alJl201ll : 
ISemena et al.l 20141), we haye vjf = 4000 km s“'. Using the strong 
shoc k condition kTshock - 3/impV^/16, we find that kTshock ~ 19 
keV. iMukai et al.l (20031) obtained a good fit to a short Chandra 
HETG obseryation with a shock temperat ure of 20 keV, Yuasa 
et al. (2010! ) obtained keV from their fit to 

Suzaku data. iBrunschweiger et al.l (2 009) fit of Swift/BAT data 
yielded kT keV, whereas iFuiimoto & Ishidal ( 1997b 
obtained k7’shock=15.4;|;2 6 keV from the line ratios in the ASCA 
spectrum. 

We therefore use a mass of 0.79 M© and a radius o f 7x10^ cm 
for the WD. The preshock density is aboye 10*^ cm“^ dLuna et alJ 
2010) which, combined with the low magnetic field strength (~ 
4x10"^ G), means that cyclotron cooling can be neglected. 


3. Observations 

EX Hya was obseryed with Chandra using the ACIS-S/HETG 
combination for 496 ks. The obseryation was obtained in four 
segments (ObsIDs 7449; start time 2007 May 13 22:15:35 UT, 
exposure time 130.65 ks; 7452: start time 2007 May 17 03:12:38 
UT, exposure time 49.17 ks, 7450: start time 2007 May 18 
21;56:57UT, exposure time 162.73 ks; and7451: start time 2007 
May 21 14:15:08 UT, exposure time 153.07 ks). We reduced 
the data using the Chandra Interactiye Analysis of Obseryations 
(CIAO y.4.1) software package. Beyo nd the reduc tion of the 
Chandra data set already described in iLuna et al.l (2010) we 
combined the +1 and -1 orders of the HEG and MEG arms us¬ 
ing the script add_grating_orders. The HEG+1 and MEG+1 
spectra were thus fit simultaneously. Line fluxes were measured 
by using Gaussians to represent the emission lines and a first- 
order polynomial to represent the nearby continuum (see Figure 
[TJ. We estimate that the measured fluxes differ by less than 2% 


using the more complex yS-profile, which indeed was loaded with 
the RMF during the fit. 

Table [T] lists the measured fluxes of indiyidual H-like and 
He-like lines. For the analysis discussed below, H-like Lya and 
the He-like resonance (r), intercombination (/), and forbidden 
(/) lines are required. For O, Ne, and Mg, the / lines are not 
obseryed. The r and i lines of Ar cannot be resolyed by the 
HEG and these two lines are measured together using a single 
Gaussian component. Because the emissiyity of the / line of Ar 
XVII 43.994 peaks at logT -7.3, while the emissiyities of S 
XVI 43.991 and 43.992 peak at log T -7.4 (AtomDB), the mea¬ 
sured flux of the line at 43.99 consists of the summed contribu¬ 
tion of these three lines, and the models take this into account. 
Since the He-like r, /, and / lines of Ca XIX cannot be resolyed, 
we used a single Gaussian to measure these lines. Finally, the i 
and / lines of Fe XXIV cannot be resolyed and thus we mea¬ 
sured thes e two lines togeth er using a single Gaussian. As dis¬ 
cussed by lLuna et al.l (20IQ) the broad wings in the residuals to 
the fits of some of the emission lines, in particular O VIII, are 
due to emission by the photoionized preshock accretion flow. 
Dielectronic recombination transitions, as well as other satellite 
lines located close enough to the He and H-like lines that cannot 
be resolyed at the HEG spectral resolution were included in the 
measured fluxes and in the models described below. We found 
that they contribute to a small fraction ( a few percent) of the 
line fluxes. 


4. The postshock structure 

4.1. Determining the shock temperature 

To ayoid the complications of modeling the complex absorber 
in an X-ray global fit, we exploited the high quality of the 
data set and filtered the data to obtain a spectrum at spin-phase 
maximum when the obseryer has a yirtually unblocked yiew 
of the accretion flow. We diyided the data into fiye spin-phase 
bins of width A(p(,j-0.2 refers to the phases at the spin 
period of 67.03 min) and fit the cooling flow model ayailable 
in XSPEC (mkcflow) to the fiye spectra. We had already fil¬ 
tered in the binary phase range <^ 98 = 0 . 0 - 0.6 (<^98 refers to the 
phases at the orbital period of 98 .26 min) to ayoid t he absorp¬ 
tion by the accretion disk bulge (IHoogerwerf et all 2006) As 
expected from obser yations at EUV wayelengths (iHurwitz et al.l 
ll997l: lM auchel 199% and from our spin-phase filter, the result 
of fitting the data in the spin-phase range 067=0.9-1.1 yields 
a yalue for the absorption column Ah smaller than 10^*^ cm“^ 
and M - 2.18;|;gg[ x 10“" M© yr“'. This fit yields a yalue 
of A:7’shock=25.5;|;Q 5 keV. If we allow the elemental abundances 
to yary (using the XSPEC model vmcflow), we get a temper¬ 
ature k7’,hock=19.7!°;2 keV, M = 2.55+°;°} x 10-" M© yr"', 
and abund ances (in un its of the solar abundances of Anders 
<& Greyesse 19891) of 0=0.92+°;°^, Ne=1.31+° °°, Ar=0.93+° g, 
Ca=1.07+Q3.}, Ee=0.61+Q}}}; Mg, S, and Si are consistent with 
sol ar yalues (these abundan ces agree with the yalues deriyed 
by lEuiimoto 8l Ishidal 19971) Although these abundances look 
less anomalous in unit^of other modern solar abundance deter¬ 
minations (e.g;JLodders_e^_alJ_2009Y we retain the Anders & 
Greyesse (19891) scale because AtomDB and XSPEC use them 
as defaults. The temperature deriyed is thus that predicted using 
eq.[T]and the strong shock condition with a WD mass of 0.79 M©, 
in agreement with yalues determined dynamically. Therefore, 
throughout the following sections, we adopt a shock tempera¬ 
ture kT =19.7 keV and the abundances deriyed from the vmcflow 
model. 
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Fig. 1. Observed emission-line profiles with overlaid Gaussian models of (a) O VII A21.6 and O VII /12I.8 (MEG +1); (b) O VIII 
dl8.97 (HEG +1); (c) Ne IX dl3.45 and Ne IX dl3.55 (HEG +1); (d) Ne X A12.13 (HEG +1); (e) Mg XI 49.16 and Mg XI 49.23 
(HEG +1); (f) Mg XII 48.42 (HEG +1); (g) Si XIII 46.64, Si XIII 46.68, and Si XIII 46.74 (HEG +1); (h) Si XIV 46.18 (HEG+1); 
(i) S XV 45.04, S XV 45.06, and S XV 45.10 (HEG +1); (j) S XVI 44.73 (HEG +1); (k) Ar XVIII 43.73, Ar XVII 43.944 -h Ar XVII 
43.961 -I- Ar XVII 43.965, and Ar XVII 43.994 -i- S XVI 43.992 -H S XVI 43.991 (HEG +1); (1) Ca XX 43.02 and Ca XIX 43.17 
(HEG +1); and (m) Fe XXVI 41.78, Fe XXV 41.85 and Fe XXV 41.86 (HEG +1), where the HETG data and model prohles are 
shown with black and red histograms, respectively. Below each panel, we plot the ht residuals. Where overlapping, th e HEG+1 and 
MEG+1 spectra were ht simultaneously; HEG+1 is shown for illustrative purposes. As discussed in iLuna et al.l (201C0 the broad 
wings in the residuals to the hts of some of the emission lines — O VIII in particular — are due to emission by the photoionized 
preshock accretion How. 
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Wavelength (Angstrom) 

Fig. 2. Chandra MEG +1 order X-ray spectrum between spin-phases <p(,j-Q.9 and 1.1 (i.e., the range of phases where the observer 
has the most unextinguished view of the accretion column) fit with an XSPEC variable-abundance isobaric cooling flow model 
(vmcflow). The absorption column, Ah, is smaller than 10^° cm“^, as expected, while the shock temperature is k7’shock=19.7);g j 
keV. 


4.2. The isobaric cooiing fiow. Derivation of a 

physicaiiy-based emission measure distribution 

In the simplest scenario, the plasma in the accretion column ra¬ 
diates its excess energy in an isobaric flow where the cooling is 
determined by the radiative cooling function A{T) (Eig.[3]l. The 
cooling time can be expressed as 



5kT 

“ nAiT) ’ 


( 2 ) 


where the density n - hqTqIT, with the 0 indicating values just 
beneath the shock. The rate of temperature change can be ex¬ 
pressed as 

dT _ noToAjT) 

dt ~ 5kjiT ’ ^ ’ 

where k^ is the Boltzmann’s constant. 

The speed in the cooling flow v = dxjdt can be written as 


V = 


vo«o 

n 


(4) 


from mass flux conservation, where vq = Vff/4. Then, at each 
place in the column, we can write 


dx - vdt — 


Vffno SkfiT SksT^Vff 

- dT = —T- dT . 

An noToA{T) 4T^noA(T) 


(5) 



Fig. 3. Upper panel: Radiative cooling function derived using 
emissivities from AtomDB (in the low-density limit, see Section 
14.21 1 and abundances (in units of the so lar abundance of Anders 
& Grevesse 19891) of 0=0.92, Ne=1.31, Ar=0.93, Ca=1.07, 
Fe=0.61 (see Section OT) . Lower panel: Line emissivities as 
a function of temperature from AtomDB, normalized by their 
peak emissivity, for the H-like (full, black line) and He-like res- 
onance-i-intercombination (dotted, red line) transitions of O, Ne, 
Mg, Si, S, Ar, Ca, and Ee. 


The emission measure (EM) is defined as EM - J n^dV. A sim¬ 
ple geometrical model for the accretion column is a cylinder, with a radius rspot and height /zsh. In a cylinder, the area element 
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dVIdx - nP' is constant, and using eq. |5]and n - noTolT, we 
have the following expression for the EM; 

EM = jjrr^ fTioVf^kB f K{T)~^dT , ( 6 ) 

where AT is the temperature range from the shock temperature 
to the minimum temperature of the gas when it reaches 

the WD surface. In practice, the value of kT„,i„ is set by the low¬ 
est temperature for which the cooling function was calculated, 
which in our case is log T-A. In fact, the assumption that the gas 
can be described in the coronal approximation, where the ions 
are in the ground state, breaks down below log T =6, but that 
region makes no significant contribution to the X-ray emission 
lines discussed here (see Fig. B- 

The flux of a given emission line is 

^ EMD(T)dT , (7) 

where A is the elemental a bundance, a constant that m od- 
ifles the solar abundances (lAnders & Grevessel 1989D al¬ 
ready em bedded in e(T), t he emissivity of the ion from 
AtomDB (iFoster et al.l 201 2D, EMD(T) is the distribution of 
the emission measure as function of temperature (i.e., the 
EmissionMeasureDistribution)', note that A(7’) was computed 
with the set of elemental abundances found in the global spec¬ 
tral fit in Section l4Tl 

While most of the constants in equation|6]are unknown, they 
primarily affect the absolute intensities of the emission lines, 
so the ratios of H-like to He-lik e emission lines are mo re eas¬ 
ily compared with observations (iFuiimoto & Ishidal 199"^ The 
constants in equationj^ affect the EMD, and hence the line ratios 
indirectly, in that they determine the shock height and therefore 
contributions of the change in gravitational potential and com¬ 
pression due to dipole geometry of the postshock flow. In the 
present case, we neglect the dipole geometry and assume a shock 
height much less than the WD radius. From equation|2] we can 
express the H to He line ratio as 

_H _ / eLyqjT) X A(T)-^dT 

f (^He, + ^He, + ^Hcf) X A{TY^dT 

where the subindices r, i, and / refer to the resonance, inter¬ 
combination, and forbidden transitions, respectively, and Lya 
includes contributions from both components of the doublet. 
Using equations 12 and 0 with k7’shock=19.7 keV, we calculated 
the H/He line ratios for O, Ne, Mg, Si, S, Ar, Ca, and Fe and 
compared them with the ratios measured in the Chandra data 
(see Table |2]). The observed ratios of the inten sities of Fe XVII 
and Fe XXII lines dMauche et akllTOOll 2003i) indicate that the 
density in the accretion column is very high, with n > 10^"^ cm“^. 
At these densities, the sum of the flux of the i and / lines is nearly 
constant, with the / line flux suppressed in favor of the i line flux 
as the density increases. Thus, the lack of observed / lines for 
O, Ne, and Mg should have a negligible effect on the sum of 
the observed He-like line fluxes. The r line is very weakly den¬ 
sity dependent. Using the sum also reduces density-dependent 
effects. We use the public AtomDB emissivities, computed for 
cm“^, and thus although the flux of the forbidden lines at 
this low density is much higher than the flux of the intercombi¬ 
nation lines, their sum is density independent. 

In Fig. |4] we plot the measured line ratios and model pre¬ 
dictions. The isobaric cooling flow model overestimates the line 


ratios for ions with Z < 14. As the emissivities of H- and He-like 
ions peak at different temperatures (Fig. O, they are emitted at 
different heights in the accretion column, and thus the observed 
disagreement in the line ratios suggests that the temperature and 
density structures of the accretion column are not adequately 
described by this simple model. Although the measured ratios 
were not corrected for absorption, the difference in the absorp¬ 
tion cross section between the H- and He-like lines of the same 
ion is so small that low-Z lines, such as O VII421.60 and O VIII 
418.97, would have to be attenuated by orders of magnitude to 
change the line ratios by factors of 2. 

To determine which physical mechanism could be missing 
in our simple model, we need to determine whether the model 
fails to reproduce the H-like line fluxes, the He-like line fluxes, 
or both. We can use equation [T] to calculate the line fluxes if 
we know the values of the constants in equation|2 the radius of 
the accretion spot rspot and the free-fall speed Vff, which in turn 
is determined by the shock temperature. We assume Vff as de¬ 
rived in Section l4~n and rspot = -\/MInvffnomu. The best match 
between data and model is obtained assuming an accretion rate 
M=1.74xl0“^' Mg yr“' (note that Beuermann et al. 2003 es¬ 
timated M=2.94xl0“" Mg yr“' for Mwd = 0.79 Mg). The 
cooling flow model closely reproduces the H-like lines fluxes, 
whereas it underestimates the He-like lines fluxes (Table[T]|. 

Using the specific accretion rate of a-3 g cm^^ s~^ and 
/=1.6xl0“^, as recently obtained by ISemena et al.l (2014D, re¬ 
sulted in both H- and He-like fluxes overestimat ed by more than 
an order of magnitude. IFuiimoto & Ishid^ t 199711 found that their 
cooling flow model fit the line ratios of high-Z ions (Z > 14), 
assuming that the bottom of the accretion column has a tem¬ 
perature kTB-Q.65 keV, preventing any cooling below this tem¬ 
perature. In this scenario, the H-like line intensities are well re¬ 
produced but the He-like lines of low-Z ions (Mg, Ne, O) are 
significantly underpredicted. 


4.3. Magnetic dipole geometry 

The simple cooling flow model (Model A in Table B makes a 
number of approximations, neglecting photoionization, depar¬ 
tures from constant pressure, departures from planar geome¬ 
try, variations of the gravitational potential with height, resonant 
scattering, thermal conduction, electron-ion temperature differ¬ 
ences, and Compton cooling (see Appendix A). The formulation 
of Model A cannot be easily modified if we want to include other 
mechanisms. We th erefore developed a modified version of the 
model presented by ISuleimanov et’aP (2()()5D using the radia¬ 
tive cooling function derived from AtomDB emissivities with 
the abundances derived i n Section 14. 1 1 an d a magnetic dipole 
geometry as described by ICanalle et aU (2005D We explored a 
range of specific accretion rates, a, between 0.2 and 2.0 g s“' 
cm“^, resulting in a range of shock tem peratures from 15.9 to 
21.3 keV. The models usef?in = 2.7 Rwd (IRevnivtsevet al.l201 lb 
ISemena et al.ll2014ll . Tables [1] and |2]list the resulting fluxes and 
line ratios from this model (Model B), respectively, with a shock 
temperature of 19.7 keV, which matches the shock temperature 
derived from the WD mass of 0.79 Mg, a specific accretion rate 
of 0.6 g s“' cm“^, and an accretion rate M=l.74x10“'* Mg yr“'. 
The line ratios are higher than the measurements, especially for 
low-Z elements. Line fluxes from H-like ions are greater than the 
measurements, while He-like lines are in closer agreement with 
the data. Although a model with a lower accretion rate would 
help to match the H-like line fluxes, there is a trade-off between 
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Fig. 4. Measured ratios of the H-like Lya line to the sum of the He-like resonance, intercombination, and forbidden lines from the 
Chandra spectrum of EX Hya {black diamonds with error bars). The histograms show the results from the models discussed in 
Sections [4.2l to l4.5l For the empirically-derived EMD, the error bars indicate the 90% conhdence level (Sect. 14761) 


matching these and the He-like lines fluxes and the continuum, 
which tends to be underestimated by models with very low M. 

4.4. Varying the elemental abundances in the cooling 
function 

In the analysis above, bremsstrahlung dominates the cooling be¬ 
tween logTshock and about logT = 7.3, since it gives the ob¬ 
served EMD and the relative H-like ion intensities (Table [TJ. 
Given that the low-Z He-like line fluxes are underpredicted, one 
possibility is that the cooling rate below log T — 7.3 is lower 
than the adopted cooling rate. One way of reducing the cooling 
rate below logT = 7.3 is to reduce the elemental abundances, 
thereby reducing the emission-line cooling from K-shell Mg, Si, 
and S and L- shell Fe, which d ominate the cooling function for 
logT < 7.3 (iFoster et al.l 2012h A model with reduced abun¬ 
dances can indeed reproduce the observed line ratios if the el¬ 
emental abundances are about ~10% of the solar abundances 
(see Model C in Tables [T] and |2] and Fig. 0]). However, since the 
bremsstrahlung continuum is not sensitive to these abundances, 
the line-to-continuum ratios of the H-like and He-like lines are 
drastically reduced and the observed spectrum is far from ade¬ 
quately ht. 

4.5. Cooling model with thermal conduction 

Another possibility is that extra heat is added in the regions of 
the column that are closer to the WD surface, thereby increasing 
the fluxes of the low-Z He-like lines. We discuss several possi¬ 
ble mechanisms in Appendix A, along with several mechanisms 
that could potentially alter the line ratio diagnostics. We con¬ 
clude that photoionization, steady-flow departures from pressure 
equilibrium, resonant scattering, Compton cooling, and electron- 


ion temperature differences are unlikely to explain the discrep¬ 
ancy. A dipole geometry and the gravitational potential differ¬ 
ence between the top and bottom of the accretion column also 
fail to explain the discrepancy, as shown by the models discussed 
above. An additional possib i lity is the rmally unstable cooling 
(llmamura & Chevalieiill984l : llnnesl 19921) Thermal conduction 
seems promising in that it transfers heat from the region above 
10^ K, where the H-like lines are formed, to the cooler regions 
where the He-like lines of O, Ne, Mg, and Si are produced. 

In clusters of galaxies, thermal energy from outer regions 
can be transported to the central cooling gas by conduction, but 
its importance i n this context has been considered doubtful by 
several authors (|B innev & Cowiel l 19811 : iTucker & Rosneil[l983i : 

[Bertschinger & Meiksinl 19861) Early studies assumed that the 

magnetic fields strongly suppress conduc tion perpendicular to 
the held. Later. iNaravan & Medvedevl (200lh showed that a tur¬ 
bulent magnetohydrodynamic medium could support a conduc¬ 
tive heat flux, which would play a signific ant role in the energy 
equation. In nonmagnetic accreting WDs. iLiu et al.l (2008h pro¬ 
posed that the X-ray spectra of dwarf nova can be fit with a 
model where the material in the innermost portion of the ac¬ 
cretion disk is evaporated into a flow, forming a corona that sur¬ 
rounds the WD. They found that only small conductivities, of 
about 1% of the Spitzer value, are able to provide acceptable 
fits, presumably because the strong shear in the boundary layer 
produces magnetic fields nearly perpendicular to the temperature 
gradient. 

In magnetic CVs, cooling flow models generally ignore 
thermal conduction because its effect is modest in the high- 
temperature part of the flow, reducing th e emission measure by 
a few percent above a few times 10^ K (llmamura & ChevalieJ 
19841) However, that energy is comparable to the thermal energy 
of the gas below 10^ K. The combination of lower temperature 
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Fig. 5. Emission measure distributions (EMDs) derived from models of the cooling of the postshock plasma discussed in Sections 
l4.2l to l43] Error bars indicate the 90% confidence level of the empirically-derived EMD (Sect. l4~6]) . 


and higher radiative cooling coefficient creates steep temperature 
gradients in that region, and the divergence of thermal conduc¬ 
tion is more important. If a large fraction of the conduction flux 
were deposited between about 2x10® and 10^ K, it would al¬ 
leviate the discrepancy in the H-like to He-like ratios of low-Z 
elements. 

We computed simple models of a constant pressure flow 
cooling from a shock temperature of 19 keV with a constant 
mass flux. That is, n is inversely proportional to T, and v is in¬ 
versely proportional to n (eq. |4]l. The choice of density is not 
important because the cooling rate is proportional to the density 
squared, while the cooling length scale is inversely proportional 
to density, and the divergence of the conductive flux is propor¬ 
tional to the inverse of the distance scale squared. We describe 
the previous discussion in the following equations; 

Pmi «e > (9) 

dx oc , (10) 

VFc dx~^ cc nl, ( 11 ) 

where Fc is the standard Spitzer conductive flux and P^ad is the 
radiative cooling rate. Thus, we have a one-parameter family of 
models provided that the shock jump conditions determine the 
relation between flow speed and maximum temperature. 

We solved the equation for thermal energy by specifying a 
temperature grid and determining the distance needed to cool 
from one temperature to the next. If thermal conduction is ig¬ 
nored, we obtain a simple linear equation (see Section l4^ . With 
thermal conduction, the equation is cubic. 

If Af, is the time needed to cool AT/ = T, - r,+i and Vi is the 
flow speed, then the thickness of that cell is 

( 12 ) 


and the contribution of that cell to the emission is n^Ax, . We find 
Af/ = (StiikBATiAxi + AFci)l(nfAAxi), (13) 

where AFc, is the difference in thermal conduction fluxes enter¬ 
ing and leaving cell i. That yields the equation 

nfAAx] - StiiVkBATiAxj - VEc/Ax; -i- VKTf^ATt = 0, (14) 

where k is the thermal conduction coefficient. 

The effect of conduction is to transfer a few percent of the 
energy from the hot part of the flow to the cooler regions, in¬ 
creasing the emission measure at lower temperatures. Nearly all 
the conductive flux is deposited below 10® K, however, and the 
EMD at higher temperatures is almost indistinguishable from the 
model without conduction. 

4.6. The empirical EMD 

We can empirically derive an upper limit on the EMD us¬ 
ing the observed emission line fluxes, their tabulated emissivi- 
ties, and a fitting algorithm. To derive the empirical EMD, we 
use a Markov-Ch ain Monte Carlo proces s as implemented in 
PINTofALE (PoA: li^shvap & Drakd 20001 A detailed explana- 
ti on of the method, its poten tial, and its limitations can be found 
in lKashvan & Dra^ ( 19981) Given the nature of the method, the 
shape of the EMD depends on the emission lines chosen for fit¬ 
ting. Basically, the EMD is obtained by fitting the minimum of 
the loci curves. The loci curves are constructed by measuring the 
fluxes of emission lines using the relationship 

EM(T) = , (15) 

AediT) 

where d is the distance to the source, is the flux of the emission 
line in ergs s * cm“^, A is the element abundance, and Q(r) is 


Ax/ = V/Af/, 
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the line emissivity as a function of temperature in units of ergs 
cm^ s ^ Although the derived EMD depends on the assumed 
elemental abundances, PoA has the option of reconstructing an 
abundance-independent EMD using the measured line ratios of 
H- and He-like lines of O, Ne, Mg, Si, S, Ar, Ca, and Ee. The 
resulting best-ht solution for the EMD is shown in Eig. 13 while 
the range of fluxes and line ratios are listed in Tables [T] and |2] as 
Model D. The overall trend of the empirical EMD agrees with 
the physically-based EMDs derived above. The empirical EMD 
mostly differs from theoretical EMDs in temperatures where no 
emission lines were measured (e.g., log T > 8.1 or log T < 6.6) 
and there are only a few continuum, line-free regions (e.g., 2.2- 

3.2A). 


5. Discussion and conciusions 

We used the high-quality spectra obtained with Chandra HETG 
from the intermediate polar EX Hya to test models for the cool¬ 
ing of the postshock flow, often applied to the X-ray emission 
from the accretion column of magnetic CVs. We have found 
that the fluxes of the H-like lines are well reproduced by mod¬ 
els of an isobaric cooling flow, while fluxes of low Z He-like 
lines are underestimated, suggesting that either the cooling in 
the low-temperature region of the accretion column does not re¬ 
move heat efficiently, or there is an additional source of heat. On 
the other hand, the X-ray data from clusters of galaxies show less 
low-temperature plasma than expected from the similar isobaric 
cooling flow scenario. Thus in galaxy clusters there must be an 
extra heat source that is strong enough to prevent the plasma 
from cooling below 10^ K (e.g.. lZhuravleva et al.l 20141) 

We studied different variables that could affect the shape of 
the EMD and the emission fine fluxes, such as a dipole geometry, 
gravity, and chemical abundances. Neither a dipole geometry nor 
the gravitational potential difference between the top and bottom 
of the column are able to modify the EMD in such a way as to 
match the observed fine fluxes and continuum simultaneously, 
and because most of the continuum is due to the bremsstrahlung 
emission at the shock temperature, it is possible that the accre¬ 
tion column is not tall enough for these effects to be significant. 
A subsolar set of elemental abundances allows us to reproduce 
the observed fine ratios, but the fine fluxes are greatly underesti¬ 
mated as are the line-to-continuum ratios. Thermal conduction, 
which has been considered in clusters of galaxies, does not re¬ 
solve the discrepancy. 

A comparison of the predicted spectra from the different 
models discussed in Section |4] and the Chandra observation is 
shown in Eig.|3 The spectra from Models A and C are not shown 
in this figure because Model A reproduces the results already 
shown in Figure|3 while Model C greatly overestimates the con¬ 
tinuum and underestimates the fine fluxes. Spectrum from Model 
B computed for an accretion rate of M=l.74x10^" Mq yr“' 
overestimates tfie continuum (Eig. 0 upper panel), indicating 
tfiat the accretion rate from this model should be lower that the 
accretion rate found in the simple isobaric cooling flow model 
(model A). Model D (Fig. flower panel) makes a reasonable 
prediction of the observed spectrum and fine ratios, although it 
underpredicts the fluxes of emission fines in the »5-10 A region. 

The appendix discusses several other possible mechanisms 
that could modify the EMD and shows that they are unlikely 
to account for the discrepancies between the model and the ob¬ 
servation for the flux of the He-like fines. Thermal conduction 
does not alleviate the problem. A possible explanation is the 
expected thermal instability of the cooling flow (Imamura & 




Fig. 6. Comparison between the Chandra MEG+1 spectrum 
(black full fine) and spectra obtained from the different mod¬ 
els discussed in Section |4] Top: Accretion column with a mag¬ 
netic dipole geometry and change of gravitational potential with 
height (blue full fine, Model B, see Section 14.31 for details). 
Bottom: EMD derived from the fit of the loci curves obtained 
from the flux measurements (orange full fine. Model D, see 
4.61 for details). Tfie model-data panels show that the 
a emission from model B overestimates the data. Given 
5 spectral models were derived assuming the same ac- 
e, M=l.74x10'^ g s“*, model B requires lower accre¬ 
tion ratels to match the data. 


Section 
continuiki 
that all t h 
cretion i a 


Chevafier ll984l : II7anger et aH 19821) so that the shock height and 
shock speed vary with time. Hydrodynamical simulations will be 
needed to determine whether such instabilities can change the H- 
like to He-like intensity ratios by the observed amount, but they 
should change the ratios in the correct sense. First, the nonlin¬ 
ear thermal instability causes the hot part of the flow to cool 
at constant density rather than constant pressure. That means 
that the energy each particle must lose is 3/2 kT rather than 5/2 
kT. Second, the energy that is not converted to radiation in the 
constant density cooling produces secondary shock waves that 
reheat the gas. The secondary shocks should have speeds of a 
modest fraction of the intial shock speed, and they could plau¬ 
sibly produce temperatures of a few million K, leading to extra 
emission in He-like fines of the lower-Z elements. This simple 
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picture is probably complicated by random variations in the den¬ 
sity of the gas striking the shock. 

Model D shows that the empirical EMD can be rather com¬ 
plex and better theoretical models will be needed to understand 
the high-resolution X-ray spectra of magnetic accreting WDs 
that will be obtained with ground-breaking future instruments 
such as the calorimeter onboard ASTRO-H. 
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Appendix A: Modifications to the cooiing modei 

There are several physical processes that must operate at some 
level in the shock-heated plasma and that would in principle af¬ 
fect the relative intensities of the H-like emission lines. 

Photoionization: X-rays from the hotter parts of the flow heat 
and ionize the cooler parts. X-ray emission from the photoion- 
ized accretion column abow the shock was identified bv Luna 
et al. (20101) in several lines in the form of very large velocity 
width components. In the shocked gas, the importance of pho¬ 
toionization is greatly reduced because of much higher density. 
We have estimated the effect of photoionization on the cooling 
fl ow by computing m odels of steady-flow shocks with the model 
of I Raymond! ( 1979D The atomic rates in this model have been 
updated over the years, but they are general ly not as accurate as 
the AtomDB rates (see lSmith et al.l 200 IL for a comparison) 
However, the percentage changes due to including processes 
such as photoionization should be reliable. Simple 3000 km s ' 
shock wave models with and without photoionization showed no 
difference in the H-like to He-like line ratios. 

Departure from constant pressure: The shocked plasma 
slows from 1/4 the shock speed to zero as it settles to the WD 
surface. That implies a pressure gradient that can influence the 
emission measure distribution. However, models with magnetic 
support cool at nearly constant density down to 10® K, and there¬ 
fore nearly constant velocity. These models show changes in the 
line ratios at only a percent level compared to models without 
magnetic support. The perpendicular component of the magnetic 
held alone can support the cooling plasma and the held is en¬ 
tirely parallel to the flow in the standard shock model. However, 
observations of supernova remnants show strong amplification 
of the magnetic held, qu ite likely by an instabil i ty that pro duces 
perpendicular field te.g.. lVink & Lamindl2003HBelll 2004), and 
this kind of mechanism may operate in accretion shocks as well. 

Thermal instability: Larger changes in the LMD could re¬ 
sult fro m the thermally unstable nature of the postshoc k cooling 
region (llmamura & Duris'^ll983l : lLanger et al.l 1982D The in¬ 
stability causes oscillations in the shock height and shock speed 
and the cooling gas can deve lop large low-pressure regions and 
subshocks (e.g., InnesI 19921) It is not obvious whether thermal 
instability would increase or decrease the He-like to H-like ra¬ 
tio, or how big the effect would be. Hydrodynamic models are 
needed. 

Resonant scattering: Terada et al. (2001, 2004) investigated 
the effects of scattering in the resonance lines in magnetic 
CVs. They found that resonance scattering produced asymmetric 
emission and modul ation on the spin pe riod in polars, but not in 
intermediate polars. iLuiimoto & Ishidal ( 19971) also showed that 
resonant scattering was not important in LX Hya. 

Electron-ion temperature difference: If the shock wave heats 
ions, but not electrons, as is observed to occur in co llisionless 
shocks in supernova remnants (iGhavamian et alJ 2013l) the elec¬ 
trons are gradually heated by Coulomb collisions. The result is 
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that the peak electron temperature is a bit lower than the nomi¬ 
nal shock temperature, and the electrons stay near that peak tem¬ 
perature longer than in a single temperature flow. Based on the 
shock models mentioned above, this increases the Fe XXVI line 
flux relative to the others, but not very strongly. 

Compton Cooling: Compton cooling can compete 
with bremsstrahlung if the density is low, reducing the 
bremsstrahlung emission at th e highest temperature s. For a 
WD temperature of 25,000 K (lEisenbart et alJ 20021) and the 
densities indicated bv the Fe XVII and Fe^Qffl lines (Mauche 
et al. 1200 L 20031), however, the Compton cooling rate is less 
than 1% of the bremsstrahlung cooling rate. Compton scattering 
is an additional cooling mechanism, which can potentially 
decrease the plasma temperature in the upper hot layers and, 
therefore, decrease the fluxes of H-like lines. However, this 
cooling mechanism is effective only for massive WDs (> 1 
Mq) and high local ma ss accretion rates (a > 1 g s“' cm“^) 
(ISuleimanov et al.1 20081) 
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Table 1. H- and He-like line fluxes (in 10 ^ photons s ' 


cm measured and derived from the models discussed in Section|4] 


Ion 

Data 

Model A“ 

Model B* 

Model C" 

Model D** 

0 VII 421.8. 

1 

^•^^-0.14 

0.92 

0.85 

0.45 

0.40-1.53 

O VII 421.6. 

^•^^-0.17 

1.11 

1.03 

0.53 

0.48-1.59 

O VIII 418.97. 

o 4Q+0.31 

8.57 

9.35 

2.24 

3.45-9.33 

Ne IX 413.55 + Ee XIX 413.551 -H Ee XX 413.535 -H 413.533 

0 

'^•'^^-0.07 

0.54 

0.53 

0.14 

0.25-0.71 

Ne IX 413.45+ Fe XIX 413.462 . 

1 47+0.08 
-0.08 

0.73 

0.75 

0.20 

0.36-1.00 

Ne X 412.13-nFe XVII 412.124 . 

4 07+^ 0^ 
^•^^-0.17 

3.36 

3.82 

0.58 

1.49-4.01 

Mg XI 49.23 + Fe XXII 49.231 . 

0 34+0 02 
'^•'^^-0.03 

0.26 

0.28 

0.05 

0.13-0.36 

Mg XI 49.16 . 

0 57+0.03 
'-0.03 

0.36 

0.39 

0.07 

0.17-0.47 

Mg XII 48.42 . 

1 51+0-03 
^•‘^^-0.04 

1.70 

1.97 

0.21 

0.65-1.72 

Si XIII 46.74 . 

^•^^-0.02 





Si XIII 46.68 + 46.663 -t 46.664 . 

0.45!“;® 

0.39 

0.40 

0.06 

0.16-0.40 

Si XIII 46.64 . 

0.83!“;® 

0.59 

0.63 

0.09 

0.24-0.59 

Si XIV 46.18 . 

1 q^+ 0.03 

2.14 

2.54 

0.23 

0.83-2.22 

S XV 45.10 +45.098 -H 45.103 . 

'^•'^^-0.02 





S XV 45.06 + 45.050 . 

n 97+0.06 
- 0.06 

0.33 

0.34 

0.03 

0.10-0.28 

S XV 45.04 . 

0.48!“;® 

0.40 

0.44 

0.05 

0.17-0.42 

S XVI 44.73 . 

1-01!!!;^^ 

1.14 

1.41 

0.11 

0.43-1.19 

Ar XVII 43.944 + 43.961 + 43.965 . 

0 14+0.03 
^•^^-0.04 

0.11 

0.18 

0.01 

0.05-0.12 

Ar XVII 43.994 + S XVI 43.991 + 43.992 . 

A 1 O+0.03 
'^■^°-0.03 





Ar XVIII 43.73 . 

A 90+0.04 

'^•^■^- 0.03 

0.19 

0.25 

0.02 

0.08-0.22 

CaXIX 43.184 + 43.211 . 

0.17!“;® 

0.13 

0.16 

0.02 

0.05-0.14 

Ca XX 43.02 . 

0.18!“;® 

0.14 

0.19 

0.02 

0.04-0.14 

Fe XXV 41.85 . 

1 96 + 0-29 
^'^^- 0.40 

1.05 

1.37 

0.21 

0.35-0.97 

Fe XXV 41.865 + 41.868 . 

1 47 + 0.42 
- 0.29 

0.68 

0.88 

0.13 

0.25-0.67 

Fe XXVI 41.78 . 

1 97+0.11 
^•^'-0,10 

0.87 

1.26 

0.16 

0.13-1.02 


Notes. Cooling flow model. Geometry of the accretion column: cylinder. No gravity. Solar abundances. Fractional area /=10“‘*, accretion rate 
M= 1.74x10“*' Mq yr“*. **’* Suleimanov et al. (2005) model. Geometry of the accretion column: magnetic dipole. Solar abundances. Accretion rate 
M= 1.74x10“*' Mq yr“*, specific accretion rate fl=0.6 g cm“^ s“*. Cooling flow model. Geometry of the accretion column: cylinder. Abundances 
0.1 solar. Fractional area /=10“‘*, accretion rate M=1.74xl0“** Mq yr“*. *'** EMD derived empirically from the measured flux ratios of H- and 
He-like lines of O, Ne, Mg, Si, S, Ar, Ca, and Fe. Error bars were determined at the 90% confidence level. 

Table 2. H/He line ratios, measured and derived from the models discussed in Section|4] 


0 

Ne 

Mg 

Si 

S 

Ar 

Ca 

Fe 


Data 


2.95±0.25 2.01±0.10 1.67±0.08 1.26±0.04 1.23±0.11 0.74±0.14 1.02±0.31 0.37±0.06 


Model A“ . 4.22 2.64 2.71 

Model B* . 4.96 2.95 2.91 

Model C" . 2.30 1.67 1.65 

Model O'* . 3.07+“;® 2.22+“;“ 1-98!“;® 


2.19 

1.68 

1.16 

1.05 

0.50 

2.43 

1.93 

1.26 

1.17 

0.56 

1.47 

1.27 

1.00 

0.94 

0.48 

1 1 +0.10 
^-0,01 

1 iro+0.13 
I.D6-021 

1 1 7+0.04 
^ ^ ' -0.09 

^•^^-0.16 

'-'■^^+0.08 


Notes. Same as Table[T] 
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